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Abstract 

We deal with an inverse problem arising in corrosion detection. The 
presence of corrosion damage is modeled by a nonlinear boundary con- 
dition on the inaccessible portion of the metal specimen. We propose 
a method for the approximate reconstruction of such a nonlinearity. A 
crucial step of this procedure, which encapsulates the major cause of the 
ill-posedness of the problem, consists of the solution of a Cauchy problem 
for an elliptic equation. For this purpose we propose an SVD approach. 
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1 Introduction 

In this paper we deal with an inverse problem originating from corrosion de- 
tection. The corresponding direct problem, which models the electrochemical 
phenomenon of corrosion, is given as follows 



(1.1) 



where Ti and I^ are two open, disjoint portions of dfl such that To = 60 \(riUr2). 
We recall that, given g £ L 2 (T2), a weak solution of problem (1.1) is a function 
ii e i? 1 (0), such that u\-p D = in the trace sense and which satisfies 
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f Vu-V? ? = / grj + f 
Jn Jr 2 JT! 



f(u)v , (1-2) 



for every r\ G i/ 1 (0) such that rj\r D = 0. 

In the above boundary value problem Q represents the electrostatic conductor, 
Ti represents the part of the boundary subject to corrosion, T2 represents the 
portion of the boundary accessible to direct inspection and is a portion of 
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the boundary where the electrostatic potential u is grounded. Such type of 
model with a specific choice of the nonlinear boundary term on Ti has been 
introduced and discussed by M.Vogelius and others [4], [11], [20]. 
The inverse problem that we want to address here is the determination of the 
nonlinearity / = f(u), when one non trivial pair of Cauchy data u\r 2 , f^|r 2 is 
available. That is. we assume that we can measure, on the accessible boundary 
r 2 , the voltage and the current density. Note that also such an inverse problem 
must be formulated in a weak sense and the Cauchy data ip, g must be taken in 
the appropriate trace spaces, see Section 4 for details. In a previous paper [2], 
we have considered the issue of stability. The main results obtained there are 
summarized in the following Section 3. Let us just mention here, that we have 
considered classes of unknown nonlinearities for which also the direct problem 
(1.1) might not be well-posed, and also that for such an inverse problem also 
the domain within R where / can be identified is part of the unknowns. 
In this paper we intend to initiate the study of a procedure for the approximate 
identification of the nonlinearity from approximate measurements of the data 
u\r 2 , f^|r 2 - K seems necessary, as a first step of such a procedure, to solve the 
Cauchy problem for u with Cauchy data on T 2 and determine the corresponding 
Cauchy data for u on the inaccessible part of the boundary Ti, where the cor- 
rosion takes place. The approximate solution of a Cauchy problem for elliptic 
equations has already been studied by many authors, just to mention some of 
the most recent contributions, [3], [6], [8], [9], [10], [13], [14], [15], [17], [18]. 
In this paper we wish to propose an approach based on the reduction of the 
Cauchy problem to the regularized inversion of a suitable compact operator by 
the use of its singular value decomposition (SVD). In Section 4 we discuss this 
approach treating the problem in the wider generality of a Cauchy problems for 
variable coefficients elliptic equations. In Section 5 we specialize the same ap- 
proach to the Laplace equation in a domain with a cylindrical geometry which 
might be well-suited to a reference conductor specimen, and to the model of 
electrochemical corrosion. We conclude in Section 6 by outlining the various 
steps of the approximate identification of the nonlinearity /. 

2 Main assumptions 

A priori information on the domain 

We shall assume throughout that fl is a bounded, connected domain in R™, 
n ^ 2 such that diam(il) D and with Lipschitz boundary dQ with constants 
ro, M. Moreover, we assume that the portions of the boundary I\ are contained 
respectively into surfaces Si, i = 1,2 which are C 1,a smooth with constants 



We also suppose that the boundary of Fi, within Si, is of C 1,a class with con- 
stants r , M. 

We introduce some notation that we shall use in the sequel, for every p > and 
i = 1,2, we set 



r ,M. 



Ul = {xett: dist(x, dn \ T 4 ) > p} 



(2.3) 




(2.4) 
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In some places, it will be necessary to isolate one privileged coordinate direction, 
to this purpose, we shall use the following notation for a point x € R™, 
x = (x',x n ), with x' e M n_1 , x n G R. 
A priori bound on the energy 

We assume the following bound on the measured electrostatic 

/ |Vu(x)| 2 s^E 2 . 
Jn 

A priori information on the boundary data 
The current flux g is a prescribed function such that 

IMIc".«(r a ) «S G , 

HslU-(iW > m > . (2.6) 
A priori information on the nonlinear term 

We assume that the function / belongs to C 0,1 (IR,1R) and, in particular, 

/(0) = and |/(u) - f(v)\ < L\u-v\ for every «,»£!, (2.7) 
A priori information on the conductivity 

In Section 4 below, we shall consider Cauchy problems for solutions to vari- 
able coefficients elliptic equations of the form div(crVu) = in fi. We shall 
assume that the conductivity tensor o~(x) = (fjj(a;))™j_ 1 satisfies the ellipticity 
condition 

n 

A" 1 < 0fc(aO&& A|£| 2 , for all x E Q and £ e R", (2.8) 

and the Lipschitz condition 

\aij{x) — <Jij(y)\ ^ K\x — y\, for all i, j = 1, . . . ,n and x, y G 0, (2-9) 
where if > 0, A ^ 1 are prescribed constants. 

From now on wc shall refer to the a priori data as to the set of quantities 
r , M, a, L, G, E, D, m, A, K. 

In the sequel we shall denote by rj(t) and ui(t), two positive increasing functions 
defined on (0, +oo), that satisfy 

for every < t ^ G , (2.10) 

u(t) ^ C |log(t)p w , for every < t < 1 , (2.11) 

where c>0,C>0,7>1,0<#<1 are constants depending on the a priori 
data only. 



potential 

(2.5) 



r](t) ^ exp 
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3 Stability result 



We review here the stability result obtained in [2]. As preliminary step, we 
evaluate the amplitude of the range of u on T\ and we prove a lower bound on 
the oscillation of u on I\. 

Theorem 3.1 (Lower bound for the oscillation). Let Q,g satisfying the 
a priori assumptions. Let u be a weak solution of (1.1) satisfying the a priori 
bound (2.5) then 

oscu ^ ??(]MU«>(r 2 , 2ro )) 

where r\ satisfies (2.10). 

Proof. [Sketch]. First we note that the (unknown) Cauchy data u\r 1} f^jlri 
can be dominated by the oscillation of MonTi. Next, we use a stability result 
for the Cauchy problem with data on Ti and obtain a bound on §77 |r 2 • We refer 
to [2, Section 3]. □ 

The main result in [2] is the following. 

Theorem 3.2 (Stability for the nonlinear term /). Letu.i £ H 1 ^), i = 1,2 
be weak solutions of the problem (1.1), with f — fi and g — gi respectively and 
such that (2.5) holds for each m. Let us also assume that, for some positive 
number m, the following holds 

IHU~(r 2|2ro ) >m>0 . (3.12) 

Moreover, let ipi = itjL , i = 1,2. There exists e$ > only depending on the a 
priori data and on m such that, if, for some e, < e < Eq, we have 

IIV'i -ip2\\L3(r 3 ) ^£ , 
ll.9i -52|U=(r 2 ) < £ , 

then 

ll/i - /2|U°°(y) < , 

where 

V = {a, 0) C [-CE,CE] , 

is such that 

r){m) 
P-a> — 

and r/,uj satisfy (2.10), (2.11) respectively. 

Proof. [Sketch]. In this case we use first a stability estimate for a Cauchy 
problem with data on T 2 . Next we show that on a suitable curve on r 1; the 
restriction of U\, u 2 are strictly monotone and their ranges agree on a sufficiently 
large interval V. On such an interval we are able to bound f\ — f% in terms of 
(wi - tt2)|r! and (^ - %?)| ri - See [2, Section 4] for details. □ 
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4 Solving the Cauchy problem 



We consider here a Cauchy problem for an elliptic equation with variable coef- 
ficients 

div(crV?i) = in n, 

u = ip on E, (4.13) 

crVit • v = g on E . 

Here the conductivity tensor cr(x) = (crij(x))"j =1 satisfies the ellipticity condi- 
tion (2.8) and the Lipschitz condition (2.9). The domain fl satisfies the same 
assumptions stated previously in Section 2 and E is an open connected portion 
of d£l which is C 1,a smooth with constant 7*0, M as it was previously stated for 
the portions Ti,^. We also suppose that the boundary of E within dft is of 
Lipschitz class with constants ro, M. Moreover, we denote with and E p the 

analogous of the sets defined in (2.3) and (2.4) respectively, with r\ replaced by 

1 I 

E. We introduce the trace spaces H 3 (E), i? Q 2 (E) as the interpolation spaces 
[ J ff 1 (E),L 2 (E)]i , [i? 1 (E),L 2 (E)]i respectively, see [16, Chap. 1] for details. We 

1 1 

shall denote the corresponding dual spaces by H 2 (E)*, i/ 2 (E)*, respectively. 
We recall that there exists a linear extension operator 

E : H*(£) -> H?(dCi) , such that E(tp) = ip on E and 
IW)H ff l(an) <G|H| ffi(E) for every ^G^(E), (4.14) 

where C > is a constant depending on the a priori data only, see for instance 
[1, Lemma 7.45]. Also we recall that the operator Eq of continuation to zero 
outside E, 

is bounded from _ff 2 (E) into 2 (dfl). Note that, by such an extension, i/ 2 (E) 
can be identified with the closed subspace of flJ (90) of functions supported in 
E C d£l. More precisely, denoting by T = <9f2 \ E and 

#o(ft,r) = {u G fZ" x (n) : u| r = in the trace sense} (4.16) 
1 

we can identify i/ 2 (E) with the trace space of Hq(£1,T) on dtt. See [16, Chap. 

1] and also, for more details, [19]. 

1 1 

Given ip G H^(jy) and g G i/ 2 (E)* we shall say that u G H 1 ^) is a weak 
solution to (4.13) if u| s = ip in the trace sense and also 

/ aVu-Vrj =< g,i]\ s > (4.17) 
Jn 

1 

for every 77 G i?Q(n,r). Here < •, • > denotes the pairing between 7J 2 (E)* and 

1 

i/ 2 (E) based on the £ 2 (E) scalar product. Our first step in the solution of the 
Cauchy problem (4.13) is the reduction to the case when ip = 0. To this purpose 
we consider the weak solution W G if 1 (SI) to the well-posed Dirichlet problem 

div(aW)=0 infl, 

W = Eip on dn. [ ' 
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Setting U = u - W and G = g - a\7W ■ v\ i e # 2 (£)*, we have that U is 
a weak solution to the Cauchy problem 

div(oW) = in Q, 

U = on £, (4.19) 

crVC/ -v = G on E. 

For every /i G _ff 2 (r)* let us consider the mixed boundary value problem 

div(crVu) = in fi, 

v = on E, (4.20) 

trVw ■ v = h on r . 

Here, denoting by Hq (0, E) the space introduced in (4.16) when T replaced with 
E, a function v G Hq({1, E) is said to be a weak solution to (4.20) if 

/ aVv ■ V77 =< h, T)\ T > for every rj G H^(fl, E). (4.21) 

It is readily seen, by the Lax-Milgram Theorem, that such mixed boundary value 

problem (4.20) is well-posed. It is also evident that, finding the appropriate 
i 

h G Hq (T)* such that aVv ■ v\ i = G, would imply that v = U and provide 

-^oo(r) 

us with the solution to (4.19). We note however, that given pq > such that 

E po has nonempty interior, it would suffice to check that for some p, < p < po, 

i 

erVu ■ v — G when both functionals are restricted to _ff 2 (E p ). In fact, this is 
a consequence of the uniqueness of the solution of the Cauchy problem when 
the Cauchy data are prescribed on E p (instead than on all of E). Thus, having 
fixed p, < p < pa, the solution of the Cauchy problem (4.19) amounts to find 

h G ff o ( r )* such that -v = G on H^ (Z p ). 
We prove the following. 

Theorem 4.1. For any p, < p < po, let T p be the operator 

T p :Hi(TY - (4.22) 
h i y crVi) • 

where v G Hq(0, E) solves the mixed problem (4.20). TTie operator T p is com- 
pact. 

Proof. By the well posedncss of the mixed boundary value problem (4.20), 
the linear operator 

h ^ v 

is bounded. 

Moreover, by a standard result of regularity at the boundary, it follows that for 
every p > 0, v G C 1 '"^^) and there exists a constant C p > depending on 
the a priori data and on p only, such that 

IMIci.»(£ p ) < Cp\\v\\ H l(Q)- 



() 



Thus the operator 

D p :H l (Q) C°' a (S p ) 
i> i > crVu • 

is bounded. Finally, since the inclusion 

ip : - ff|o(Sp)* 

is compact and T p can be factored as T p = i p o D p o 5, the thesis follows. □ 

Being T p a compact operator between Hilbert spaces, we have that there exists a 
triple {&j, hj, gftj^i called singular value decomposition, such that {a? }°°_^ is a 

non increasing infinitesimal sequence of nonnegative numbers, {hj}'^_ 1 , {gj 

i i 

are orthonormal bases for 7J 2 (T)* and for i7 2 (E p )* respectively, and moreover 
it holds 

T p hj = aPgP, for every j = 1, 2, . . . (4.23) 
W - for every j = 1, 2, . . . (4.24) 

where T* denotes the adjoint operator to T p . By the regularization theory for 
the inversion of compact operators, we have that, denoting with (•, •) i the 

#00 ( s p)* 

1 

scalar product for the Hilbert space i? 2 (£ p )*, the family of operators R a , a > 

(4.25) 



5 ^ E^^fcfDj,. 



is a regularization strategy for T p , namely 

lim R a T p h = h, for every h G #| ( r )*i ( 4 - 26 ) 
(see for instance [12, Chap. 2]). Moreover, the choice 

a{e) = s 2(1 -^ (4.27) 

where 7 is a fixed number, < 7 < 1, is an admissible one, this means that if 

1 1 
given, for every e > 0, g,g £ G i? 2 (S p )* and ft G iJ 2 (r)* such that 

<? = T p /i and \\g-g e \\ 1 <e, (4.28) 
then it follows that 

lim\\R a(e) g e -h\\ 1 =0. (4.29) 

We can return now to the Cauchy problem (4.13), when tp is arbitrary in (E). 

Let us suppose that, for every e > 0, ip e G Hz(Y>), g e G iJ Q 2 (E p )*, and let 
We G _ff 1 (fi) be the weak solution of (4.18), with tp = i/j e . Let us denote by 

•R e = R a (e){ge - crVW 7 ,, • i/| E ) + crVWe • i/| r G #rio( r )*\ wh ere i? Q and a(e) 
are the regularization strategy and the regularization parameter introduced in 
(4.25) and (4.27), respectively. We propose as approximate regularized solution 
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to the problem (4.13) the function u £ € i7 1 (SI) which is a weak solution of the 
mixed boundary value problem 



div(crVu £ ) = in fi, 
u £ = ip £ on E, 

ctVu f • v = R £ on T. 



(4.30) 



In analogy to (4.20) and (4.21), we shall call weak solution of the problem (4.30), 
a function u e G H 1 ^) such that = ip e in the trace sense and such that 



/ erVu e • V77 =< R £ , rj\ r > for every 77 e Hq (0, E) 
Jo 



(4.31) 



The well-posedeness of problem (4.30) is again a consequence of the Lax-Milgram 
Theorem. The following Theorem provides a convergence results for the pro- 
cedure of regularized inversion of the Cauchy problem (4.13) that we have just 
outlined. 

1 I 

Theorem 4.2. Let tp £ H 2 (E) and g £ i? 2 (E)* and suppose that there exists 
u £ if (fi), which is a weak solution to the Cauchy problem (4.13). //, given 

e > 0, we /iaue £/ia£ Ve G ^ 2 (E) and <? s G -ff 2 (E p )* 



/7/ 1 



limu 6 |r=u|r in ff 2 (r) 



lim erVit e ■ v\ T = o~S7u ■ v\ 



#<?o(r) 



(4.32) 
(4.33) 

(4.34) 
(4.35) 



Proof. Let us observe that given S any open and connected portion of dfl, 
the following holds 

\\*VW e 'v\ s -<rVW-u\ s \\ 1 ^<(*\\W-W e \\& [a) ^C2\\E1> e -E1>\\ h 

then replacing in (4.14) ip with ip e — -0, we have by (4.32) that 

< c 3 e , (4.36) 



\o-VWe ■ v\ s - cr\7W ■ u\ s \ 



where ci, C2, C3 > are constants depending on the a priori data and on S only. 
Thus by (4.36), with S = E p , and by (4.33), we have that 



lim 



crVW e • v\ 



aVW ■ v\ 







(4.37) 



Moreover, we have that (4.35) follows by applying (4.36) with S = T, (4.29) 
with g e replaced with g £ — aVW e ■ z^|s p and (4.37). Indeed, we have 



\aVu-v\ T - aVu £ -v\ T \\ 1 < 



+ 



\aVW ■ v\ T - aVW E ■ v\ r \ 
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Finally, by a standard trace inequality 

IMr - u £ |r|| H i (T) < c 4 \\u - u E \\ H i (n) < 

< c 5 (\\*Vu ■ v\ T - aVu e ■ "Irll^j. + U ~ ^H^(S)) ( 4 ' 38 ) 

where 04,05 > are constants depending on the a priori data only, then (4.34) 
follows by recalling (4.35) and from (4.32). □ 



5 A special case 

Let D be a bounded domain in R™ -1 , with Lipschitz boundary dD with con- 
stants ro, M. From now on we shall consider this special choice of O 

fi = Dx(0,l), r 2 = D x {0} , ri = D x {1} , T D =dDx (0, 1). 

In the following we will denote by Afc, ifk, k = 1, 2, . . . , the Dirichlet eigenvalues 
and eigenfunctions of - A on D, namely 

-Aip k = \ k tp k in D, . . 

We recall that the family {<fik}%Li is an orthogonal basis in L 2 (D) and also in 

Hq(D). In the following we shall refer to the {fkW^Li as the basis normalized in 

1 

the L 2 (D) norm. We have that ip £ Hq (D) if and only if its Fourier coefficients 

= I ^k (5.40) 

JD 

satisfy 

00 

^AaM 2 <oo (5.41) 
fc=i 

1 

and that, as a norm on Hq (D) we can choose 



1 

Moreover, h S Hq Q (D)* if and only if, its Fourier coefficients 

/i fe =< /i, > , (5.43) 

satisfy 

00 

]TA fc -H fc 2 <oo (5.44) 
fc=i 
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and the norm on iJ 2 (D)* turns out to be 




1%. = I 

Here < •, • > denotes the pairing between Hq (D)* and iJ 2 (Z?) based on the 

L 2 (D) scalar product. Note also that {Xk ^fk} and {Xk^fk} constitute or- 

i i 

thonormal bases for Hq (D) and Hq (D)* respectively. 

Due to the cylindrical geometry of f2, we remark that we can identify the spaces 

Hg (Ti), H^ (Ti)*, i= 1,2, with H^(D), H^{D)* respectively. Furthermore, 

i 

as noted in Section 4, we can identify Hq Q (Ti) with the trace space on dfl of 



Hq(Q,, r) when T =(r 2 U Tjj), and the same holds when the roles of Ti and r 2 

are exchanged. 

i i 

Let xp £ i7 2 (r 2 ), g £ -ff 2 (r 2 )* and let us consider the following Cauchy 
problem with auxiliary homogeneous condition on 

Au = in fi, 

u = ip on r 2 , 

du (5.46) 
— = g on T 2 , 
ov 

u = oiiTi). 

We shall say that u is a weak solution to the problem (5.46) if u\ o = Eq(i(>) 

(r 2 ur D ) 

in the trace sense and if 

f Vit- Vf]^<9,v\r 2 > for every V £ H^(Q,(T\UT^)). 

Here Eo(ip) denotes the extension of i/j by zero outside r 2 and < •, • > denotes 

i i 

the pairing between i? 2 (r 2 )* and i/ 2 (r 2 ) based on the L 2 (r 2 ) scalar product. 
We shall use a strategy similar to the one discussed in Section 4, but with some 
slight variations, suggested by the presence of the portion of the boundary 
where u = 0. As before, we reduce the problem (5.46) to the special case when 
%l> = and introduce the well-posed Dirichlct problem 

Au = in 0, 

v = £ an.Ti, ^ (5.47) 

v = o on (r 2 ur D ), 

i 

where £ is a prescribed function in i? 2 (ri). To this purpose, in analogy with 
(4.18), we consider W £ ff 1 (0) as the weak solution to the Dirichlet problem 

AW = in 0, 

W = ip onT 2 , ( 5 _ 48 ) 

w = o on (r x ur D ) . 

The difference U = u — W shall satisfy (5.46) with -0 = and g replaced with 

G = g-m i . 

9v Wo(r 2 ) 
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Note that the well posed boundary value problem (5.47), will take the place of 
(4.20). We intend to invert the map 



dv 
av r 

in order to solve the Cauchy problem. It is 
identification of the trace spaces on T^, i = 
D. 



Lemma 5.1. Let T be the operator 



(5.49) 

convenient at this stage to recall the 
1, 2 with the corresponding ones on 



T-.HUD) - H$ (D)* (5.50) 
dv 



oi 



(5.51) 



where v is the weak solution of the problem (5.47). Then T extends to a compact 

r i i *i 00 

and self-adjoint operator on L 2 (D), such that < — Afc 2 (sinh(Afc 3 )) j <fik \ are 

its eigenvalues and eigenf unctions respectively. The singular value decomposi- 

i i 
tion of T : Hq Q (D) — * Hq Q (D)* is given by 

{-(sinh(A fc 5))- 1 ,A fc -^ fc ,A fc ^ fc }- 1 . (5.52) 

Proof. Let us first observe that the operator T is well defined since the 
problem (5.47) is well-posed. In this special setting we can represent the solution 
v of (5.47) by separation of variables, namely 

v(x',x n ) = V" — smh(\ k * x„)ipk(x') (5.53) 

k=1 sinh(A fc 2 ) 

where {£, k }kLi are the Fourier coefficients of £ with respect to the L 2 (D) basis 
{<Pk}kLi' After straightforward calculations we have that 

(5.54) 

thus the operator extends to a self-adjoint operator on L 2 (D) and since the 

eigenvalues are infinitesimal we conclude that T is compact as an operator from 

i 

L 2 (D) into L 2 (D). Moreover, since Hq (D) is continuously embedded in L 2 (D) 

and L 2 (D) is continuously embedded in H^(D)* , also T : H^(D) -» H^ {D)* 
is compact and its SVD turns out to be (5.52). □ 



As a consequence of the above Lemma 5.1, we obtain that the family of operators 
R a :H^{D)* — ► if| (.D), such that 



R a {G) = Y, (-sinh(A^))(G,^) i t ^ (5.55) 
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where [ik = (sinh(Afc 3 )) -1 , is a regularization strategy for T and the choice 
(4.27) for the parameter a is still admissible. We are in the position now to 
present the regularized approximate solution for the following special case of 
the problem (5.46). That is, given G G H^ (T 2 ), 

AU = in Q, 

u = o on r 2 , 

dU (5.56) 
— = G on T 2 , 
ov 

U = on T D . 

In this section we shall denote by [r] the integral part of the real number r. 

i 

Theorem 5.2. For every e > 0, let G £ G Hq (T2)* and suppose that there 
exists U G H 1 ^), which is a weak solution of the problem (5.56). If we have 

\\Ge-G\\ i 

then for every choice of 7,0 < 7 < 1, the function 
[log^- 1 )]"- 1 

U E (x',x n )= (-Xk~^G k , e )smh(X k ^x n )<p{x') (5.57) 

fc=i 

where {Gk, 6 }k%i are the L 2 (D) Fourier coefficients of G £ (according to the 
formula (5.43) ^, satisfies 

]imU £ \ Tl = U\ ri in fr^(Ti) . (5.58) 

Proof. Since the one defined in (5.55) is a family of regularizing operators 
and since the choice (4.27) is admissible, we have that 

limp a(e) (G e )-C/| ri |L 4, n = 0. (5.59) 

By the asymptotic bounds of the eigenvalues of the Laplace operator (see for 
instance [5, Chap. 12]) we have that there exist constants c, C > depending 
on the a priori data only, such that 

cfc^r <; \ k <; Cfc^r, fc=l,2,... . 

Thus it follows that the integer k such that ^fc > a(e) is of the order [log (e 7-1 )]" -1 . 
Moreover, since 

the thesis follows immediately by (5.59). □ 

The following Corollary 5.3 provides us with the approximate regularized solu- 
tion to the Cauchy problem (5.46). 

1 1 

Corollary 5.3. For every e > 0, let ip e G Hq Q (T2), g e G i? 2 (r2)* and suppose 
that there exists u G H l {il) which is a weak solution of the problem (5.46). If 
we have 

— -011 i < e (5.60) 
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he-q\\ i <e 



then for every choice of 7, < 7 < 1, the function 

U e {x , Xyi) 

fc=l 

sinh(A fe *(l -x n )) 



[log^- 1 )]- 1 

E ("Aft 2 Gfc. e )sinh(A fe 2 a;„)^ fe (x') + 
fc=i 



(5.61) 



(5.62) 



/c=i 



sinh(A fe 2 ) 



<Pk(x'), 



where 



Gu 



9k,,. 



i/j k:£ \k 2 coth(A fc 2 ), A; = 1,2, 



(5.63) 



{V'fcjelfeLi) {fffe,e}fc^Li are ^ e L 2 (D)- Fourier coefficients of ip e and g e respec- 
tively, is an approximate regularized solution of (5.46). Moreover, we have 



limuelr! = w|rj in #00 ( r i) 



lim 



9u 
9z> 



(5.64) 
(5.65) 



ri 



Proof. Let be the solution of (5.48) with ijj = ip s , respectively. Thus we 
can decompose u = U+W where U is the solution of (5.56) with G = g— 7^7" |r 2 - 
Moreover, by (5.60) we have 



dW e dW 
dv dv 



< -W\\ m[ a) ^C 2 \\E^ £ - E^\\ H ^ [dn) < 



< Call^-VlLi^.^C-se 



(5.66) 



where C, > 0, i = 1,2,3, are constants depending on the a priori data only. 
Thus denoting with G e = g e — ^^|r 2 7 (5.61) and (5.66) leads to 



"00 ("2)* «00( r 2j 



^00 (r 2 )* 



where C > is a constant depending on the a priori data only. By (5.59) in 
the proof of Theorem 5.2 and recalling that W = on T±, we have 



lim||i? Q(e) (G £ )- U | ri ||^ (ri) 







(5.67) 



Finally, let us consider the following Dirichlct problem 



Au e = in Q, 

u e = R a (e){G e ) on r b 







on r 2 , 
on r D , 



(5.68) 



we have that 

du e du 
dv dv 



^ooTi 



Ci\\u £ - u\\ H i {n) < 
^ C 5 (\\R a(e) (G e ) - u\ ri \ 



Us 



4>\\ i ) 
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where C4, C5 > are constants depending on the a priori data only, thus by 
(5.67) and by (5.60) 

du E du 
dv dv 



lim 



0. 



After straightforward calculations, (5.64) and (5.65) follow. 



□ 



Thus, for a given error level e > 0, the regularized solution of the Cauchy 
problem (5.46) is given by (5.62) and in particular we obtain the following 
formulas for the Cauchy data on Ti as follows 



u e\r 1 



[log^" 1 )]* 

£ 

k=l 



(A fc =Vfc, e coth(A fe 2 )-g k ,e)h 2 smh(\ k ix n )(p{x') (5.69) 



du £ 
dv 



[log^- 1 )] — 1 

(A fe -^ M coth(A fe ~^) - 5 fe;E ) cosh(A fe ?)tp(x') + 

fe=i 



E -- 



fe=l 



sinh(A fe 2 ) 



'/'ib (a;') 



(5.70) 



where the coefficients ^k,e and <?fc. e , with fc = 1,2,..., are the Fourier coefficients 
of tp E and g e , with respect to the L 2 (D) basis {v/clfc^i- 



6 Conclusion: a procedure for reconstruction 

We now return to the inverse problem of determining the nonlincarity / in (1.1) 
when the measurement u\r 2 = ip is available for a given Neumann data g. First, 
we use the methods described in Sections 4, 5, for the solution of the Cauchy 
problem. In Subsection 6.1, we outline the adaptations to the method of Section 
4 needed for our corrosion problem. In Subsection 6.2 we propose a method for 
the identification of the nonlincarity / from approximate values of u|rij §7j|rv 

6.1 Solving the Cauchy problem 

• We need to solve a Cauchy problem of the form 



< 



Au 


= 


in n, 


u = 


il> 


on r 2 


du 




on r 2 


dv 


= 9 


u = 





on T D 



(6.71) 



where u £ H 1 ^), and where in this special setting we choose ip £ Hq (T2) 

1 

and we have g € L 2 (r 2 ) C iJ 2 (r2)*. The procedure introduced in Section 

o 

4 can be applied by considering a = Id, E = T 2 , T =(ri UTd). Note 

that in this case, we have ip £ ^00(^2)- Therefore, it is convenient, in 
the formulation of the Dirichlet problem (4.18), to replace the Dirichlet 
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data E(ip) with Eo(ip). We consider W as the solution to (4.18) with such 
modified Dirichlet data, that is 



AW = in n, 
W = E {ip) onSfi. 



(6.72) 



Performing as before the decomposition u = U + W, we obtain that U is 
the solution to the following variant of the Cauchy problem (4.19) 



AU = 
U = 

au 

ov 

U = 



dW 
dv 



in Q, 

on r 2 , 
on r 2 , 

on T D 



(6.73) 



We can use the SVD decomposition described in (4.23). Note that here 
~E p = T% „ and v turns out to be the solution of the following problem 



Au = in Q, 

v = on T2, 
dv 

— = h on Ti, 
ov 

v = onTfl. 



(6.74) 



According to (4.25), we obtain a regularized inversion procedure for T p . 

• We obtain an approximate regularized solution to (6.71) by solving the 
analogue of the mixed boundary value problem (4.30), which in detail, 
takes the form 



Au e = 



u 

du, 
dv 

u E = 



in £1, 

on r 2 , 



-Rc(e)[ge- -g? | r2p )+ -gf | ri onr i' 

on T_d, 



(6.75) 



where ip e £ Hq (T2), g e £ Hq Q (T2. p )* arc the approximate Cauchy data 
and where W s £ H 1 ^) is the weak solution of (4.18), with <r(x) = Id and 
with E(ip) replaced by E (ip e ). Having solved (6.75) we can determine the 
approximate regularized values of u\r 1 , f^lri according to Theorem 4.2. 

We observe that if the conducting specimen has the special geometry introduced 
in Section 5, that is £1 = D x (0, 1), then the above described scheme simplifies 
to the formulas (5.69) and (5.70). 



6.2 Solving the algebraic equation f(u) = |^ 

We cannot expect that, for the regularized solution u e , the Neumann data 
on Ti is precisely constant on each level set of it e |ru as it should happen for 
the exact solution u to (1.1). Therefore, it is necessary to extract an approxi- 
mate expression of the nonlinearity / = f(u) when u s \r 1 and |ri m ay have 
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different level sets. We propose to obtain such approximate nonlinear term by 
minimizing the best fit functional defined as follows, 



By the Coarea formula, (see for instance [7, Chap. 3] ), we have that we can 
express F £ [f] as follows 



here, by u n -2 we denote the (n — 2)-dimensional Hausdorff measure. Thus, by 
formal differentiation it follows that 



Hence a candidate minimizcr for F £ is given by the following weighted average 
of ^flri on the level sets of u e |ri, that is 



We note the consistency of this formula in the limiting case when u e is replaced 
by the exact solution u. In-fact, in this case, the above formula leads to the 
correct values of / for every regular value t of u\r x - 
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